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Abstract 

We study here polaron (soliton) states of electrons or holes in a model describing carbon- 
type nanotubes. In the Hamiltonian of the system we take into account the electron-phonon 
interaction that arises from the deformation dependencies of both the on-site and the hopping 
interaction energies. Using an adiabatic approximation, we derive the equations for self- 
trapped electron states in zigzag nanotubes. We find the ground states of an electron in such 
a system and show that the polaron states can have different symmetries depending on the 
strength of the electron-phonon coupling. Namely, at relatively weak coupling the polarons 
possess quasi-one-dimensional (quasi-lD) properties and have an azimuthal symmetry. When 
the coupling constant exceeds some critical value, the azimuthal symmetry breaks down and 
the polaron spreads out in more than one dimension. 

We also study polarons that are formed by the electrons in the conducting band (or 
by holes in the valence band) in semiconducting carbon nanotubes. We show that their 
properties are more complex than those of quasi-lD ground state polarons. In particular, 
polarons in semiconducting carbon nanotubes possess an inner structure: being self-trapped 
along the nanotube axis they exhibit some modulations around the nanotube. 

1 Introduction 

Over the last few years much work has been done on studying physical properties of carbon nan- 
otubes [1, 2, 3], and boron nitride nanotubes [4]. The experimental studies of such nanosystems 
have revealed their peculiar properties that are important for practical applications [5]. 

Not surprisingly, carbon and boron nitride nanotubes are quite complex systems. Their ge- 
ometry is based on a deformable hexagonal lattice of atoms which is wrapped into a cylinder. 
Experimental and theoretical studies show the important role of the nanotube geometry: many 
properties of nanotubes can be modified in a controllable way by either varying the nanotube 
diameter and chirality, i.e. the way the lattice is wrapped into a cylinder, [1, 2], or by doping 
them with impurity atoms, molecules and/or compounds [6]. Theoretical studies of single wall 
carbon nanotubes (SWNT) [7, 8] have demonstrated the importance of the interaction of elec- 
trons with lattice vibrations [9, 10, 11, 12, 13, 7, 14, 15]. Note that sufficiently long SWNTs can 
be considered as one-dimensional (ID) metals or semiconductors depending on their diameter 
and chirality [1, 2]. The nanotubes possess a series of electron bands, which can be determined 
by ID energy dispersion relations for the wave vector k along the axis of the nanotube. 

In ID systems the electron-phonon coupling can lead to the formation of self-trapped soliton- 
like states (large polarons) which can move with a constant momentum [16]. In ID metals, due 
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to the Peierls instability [17], the energy gap appears at the Fermi level and the Prohlich charge- 
density wave is formed [18] instead of a soliton. Recent experiments [19, 20] have shown that 
even long channel semiconductor SWNTs may have very high mobilities at a high doping level. 
The posibility for the formation of states which spontaneously break symmetry in carbon nan- 
otubes has been discussed in [10, 11, 12]. In particular, large polarons (solitons) in nanotubes 
have recently been studied in [13, 21] where the long-wave approximation has been used for the 
states close to the Fermi level. However, such a description, equivalent to the continuum approx- 
imation, does not take into account some important aspects of the crystallographic structure of 
the system. 

In this paper, first, we consider the ground states of a quasiparticle (electron, hole or exciton) 
in the zigzag nanotube system and, second, we study the polaron states of an electron in the 
lowest unfilled (conducting) band or an extra hole (an electron deficiency) in the highest filled 
(valence) band in carbon nanotubes. For this we use the semi-empirical tight-binding model 
with nearest-neighbour hopping approximation [1]. The advantages of this method for some ID 
systems, like polyacetylene and carbon nanotubes, have been demonstrated in [22] and [1, 9, 10], 
respectively. We study a quantum system involving a hexagonal lattice of atoms and electrons 
and then perform an adiabatic approximation. Then we derive the system of discrete nonlinear 
equations, which as such, can possess localised soliton-like solutions. We perform an analytical 
study of these equations and show that, indeed, this is the case, and various polaron states can 
be formed in the system. In fact, these equations were used in [23] to determine numerically 
the conditions for the formation of such polaron states. Our analytical results on self-trapped 
states of a quasiparticle are in good agreement with the results obtained in [23] . We also study 
polarons that arc formed by the electrons in the conducting band (or by holes in the valence 
band) in semiconducting carbon nanotubes. 

The paper is organised as follows. The next section presents the model of the nanotube. The 
phonon Hamiltonian is discussed in Sect. 3, electron in Sect. 4 and in Sect. 5 we discuss the 
electron-phonon interactions. The details of the diagonalization of the electron Hamiltonian are 
presented in Appendix 1. In Section 6 we determine the adiabatic and non-adiabatic terms of the 
Hamiltonian. The corresponding zero-order adiabatic approximation then leads to the equations 
for the self-trapped electron states while the non-adiabatic term of the Hamiltonian provides a 
self-consistent test to determine the conditions of applicability of the adiabatic approximation. 
The system of equations in the zero-order adiabatic approximation in the site representation 
is derived in Appendix 2. In Sect. 7 we derive some analytical solutions for the large polaron 
ground state, and in Sect. 8 we discuss the transition to the states with broken axial symmetry. 
In Section 9 we study large polaron states in semiconducting carbon nanotubes. The paper ends 
with conclusions. 

2 Model of a Nanotube 

In this section we define the variables to describe a nanotube. Let d be the length of the side of 
the hexagons of the nanotube, R its radius and let N be the number of hexagonal cells wrapped 
around the nanotube. Then we have 

a = 27r/iV, a = dV3, b = d/2, a = 4i?sin(^), (1) 

where a is the distance between two next to nearest neighbour sites. 

To label all sites on the nanotube one can use two different schemes. They involve having 2 
or 4 lattice sites as a basic unit. The first one, used in [23], is closely connected with a unit cell 
of a graphene sheet and is based on nonorthogonal basic vectors. The corresponding labelling, 
{i,j,p), involves index i numerating the sites around the nanotube, spiral index j and index 
p = 0,1 that describes sites that have the nearest neighbours 'down' {p = 0) or 'up' (p = 1), as 
shown in Fig.l. 
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Figure 1: The two index labelling on the Hexagonal lattice. 



Note further that a hexagonal nanotube possesses two symmetries: the translation along the 
axis of the nanotube by 3d and the rotation by an angle a around the nanotube axis. Given 
this, one can use an alternative labelling scheme in which the basic unit cell is rectangular and 
contains four atoms. This scheme, also shown in Fig. 1, involves using the labelling (m,n,Q) 
where m is the axial index, n is the azimuthal index and the index g = 1, 2, 3, 4 enumerates the 
atoms in the unit cell. 

The position of any nanotube lattice site, at its equilibrium, can be described by R% given 

by 

i?° = i?(ej, sin 835 + cos 9ae) +6*2^33, (2) 

where the three-component index ae = {331,862,863} indicates the nanotube atoms, and the 
coordinates 9, being an azimuthal angle and z, being a coordinate along the tube, describe 
positions of atoms on the cyllindrical surface of the nanotube. In the first scheme 86 = {i,j,p} 
and in the second one 86 = {m,n,Q}. In a zigzag nanotube the azimuthal and longitudinal 
positions of atoms arc: 

©m,n,l = ©m,n,4 = na, &m,n,2 = ©m,n,3 = (n + ^)a; 

-1 Q — 4: 

Zm,n,e=i,2 = (3m - 1 + ^)d, Zm,n,e=3,A = (3m + IH (3) 

Although in the numerical work reported in [23], we have used the first description, the 
second one is more convenient when taking into account the boundary conditions. The azimuthal 
periodic condition /(n + A'') = /(n) is natural because going from n to n + N corresponds to a 
rotation by 27r. In the m direction, however, for a long enough nanotube, we can use the Born- 
Karman periodic conditions for electron and phonon states in a nanotube. Thus, nanotubes can 
be considered as ID systems with a complex inner structure. 

Next, we consider displacements from the equilibrium positions of the all sites of the nan- 
otube: 

R^ = Ri + U^, (4) 
where the local displacement vector can be represented as the three orthogonal local vectors: 

{/ae = Uae + Sae + V^,- (5) 



e 
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Here is tangent to the surface of the undeformed nanotube and perpendicular to the nanotube 
axis, is tangent to this surface and parallel to the nanotube axis, and is normal to the 
surface of the nanotube. Then, using Cartesian coordinates, we have 

= U^{ex cos 6ae - sin 0ae) , 
Sse = Sae(eisin6ae + e^COsGse), 

Zae = V^e^. (6) 

To write down the Hamiltonian in a compact form, it is convenient to define the formal index 
operators of lattice translations: r(), l{) and d{), which when applied to any lattice site index, 
translate the index to one of the three nearest sites. Applying these operators to the lattice 
site which has the nearest neighbour down, i.e. which in the first formulation have the index 
p = 0, they translate the index respectively to the right, left and down from that site. For the 
lattice sites which have an upper nearest neighbour, i. e. which in the first formulation have the 
index p = 1, one has to turn the lattice upside down before applying these definitions. Notice 
that the square of each of these three operators is equivalent to the identity operator. So, for 
example, moving from a lattice site to the right once and then moving to the right again, after 
flipping the lattice upside down, one returns to the starting site. In particular, we have for the 
first lattice parametrisation 

rii,j,0) = r{i,j,l) = {i,j,0), 

d{i,j,0) = {ij-l,l), d(^,i,l) = (^,j + l,0), (7) 

while for the second one, which we will use below, we have 

r(m, n, 1) = (m, n, 2), r{m, n, 2) = (m, n, 1), 

r(m, n, 3) = (m, n + 1, 4), r(m, n, 4) = (m, n — 1, 3), 

l{m, n, 1) = (m, n — 1, 2), l{m, n, 2) = {m, n + 1, 1), 

l{m,n,3) = (m, n,4), l{m,n,A) = (m,n, 3), 

d{m, n, 1) = (m — 1, n, 4), d{m, n, 2) = (m, n, 3), 

d(m, n, 3) = (m, n, 2), d(m, n, 4) = (m + 1, n, 1). (8) 

Some physical quantites, e.g. the potential energy of the lattice distortion, include central 
forces, which depend on the distance between two sites. Let us define the following lattice 
vectors connecting the atom {ae} with its three nearest neighbours (5(3e) with 6 = r,l,d for the 
right (r), left {I) and down or up (d) neighbours: 

D5^ = R5i^)-R^ = D6l + {Us^^)-U^). (9) 

When Ug. = Owe add the upper index to all quantities to indicate their values at the equilibrium 
position. Note that |-Drgg| = \Dlg^\ = \Ddg^\ = d. In the case of small displacements, i.e when 
|C/5(ae) — Usi\ d, the distance between the lattice sites is approximately given by: 

\D5^\^d + W5^, (10) 

where 

WS^ = — ^ (11) 
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are the changes of the distances between the nearest neighbours due to site displacements. The 
exphcit expressions for WSg. in the first scheme are 

Wn,j,o = ^(cos(^)(ujj,i - Uijfi) + sin(^)(si,j- 1 + Sij,o)) + ^{vi,j,i - Vij^o) 
^kjfi = (cos(^)(uij,o - + sin(^)(si_ij- 1 + Sij'o)) + - Vijfi) 

Wri,,,i = VFrij,o, Wlij,i = Wli+ij,o, Wdij,i = Wd^,,+l,o. (12) 

Because the central forces between neighbouring sites do not provide lattice stability, in 
addition to W6^, which are invariant under translations, we need also the quantities 0.6^; which 
describe relative shifts of neighbouring sites. The corresponding explicit expressions are: 

^ri,j,o = ^(cos(^)(tXij, 1 - Ui,j,o) + sin(^)(sij,i + Sjj,o)) - ^{vi,j,i - Vi,j,o), 
^kjfi = ^(cos(^)(nij-o - + sin(^)(si_ij,i + Sij,o)) - ^{vi-i,j,i - Ujj,o), 

^k,j,i = ^k+i,j,o, ^di,j,i = ^di,j+i,o- (13) 

Note, the curvature of the lattice and corresponding bond-bending in nanotubcs is an impor- 
tant factor for the lattice stability [7] and electron-phonon interaction [24] . To take into account 
this factor we choose to based our discussion on the solid angle spanned by the 3 lattice vectors 
at a given site: 

\Dr^\\DU\Dd^\ ' Id ' 
where = | sin(^) and, in the case of small displacements, 

Ci,3fi = ^ sin(|)(2t;ij,o - - t-i-ij-.i) 

- cos(-)sij_i,i + 3cos (-)si,j,o 

3 a 5 •^.ttsx. X 

+ (2 COS(-) - - cos (-))(Si_lJ,l + SiJ,l) 

+ sm(-)(-cos (-) - l)(txij,i -«i_ij,i), 

^ \/3 . ,a, , ^ . 

= — sm(-)(i;jj,o + t'i-Hij,o-2?;ij,i) 

- cos(-)sij+i,o + 3cos'^(-)sij,i 

+ (2 COS(^) - - COS (-))(Sj+ij,o + Si,jfi) 

+ sin(-)(- cos (-) - l)(«i+ij,o - Mi,i,o)- 

(15) 

It is easy to write down the corresponding expressions in the second labelling scheme. This time 
one has twice the number of the expressions as compared with the first scheme. 
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3 Phonon Hamiltonian 



We define the phonon Hamiltonian in the nearest-neighbour interaction approximation and take 
into account the potential terms responsible for the central, Vw, non-central, Vq, and the bond- 
bending, Vc forces in the harmonic approximation: 

Hph=lj:{§ + kJ2[WSl + nSl] + k,C^), (16) 

ae (5 

where M is the atom mass, k is the elasticity constant for the relative atom displacements, kc 
is a characteristic constant of the bond-bending force while is the momentum, canonically 
conjugate to the displacement Us>. 

According to the theory of lattice dynamics (see, e.g., [26]) the Hamiltonian (16) can be 
diagonalised by some unitary transformation. For the lattice labelling ae = {m, n, ^}, this 
transformation has the form 



Y^e'i>^^+--W,,r{k,i^)Qk,,,r. (17) 



Then, introducing the operators of creation, h^ ^^,^, and annihilation, bk,i/,T, of phonons 
we can rewrite the phonon Hamiltonian (16) in the standard form 

Hph = lT.{PUrPk,u,r + U;Uk,l^)Ql,,rQk,u,r) 

k,v,T 

J2 ^<^Aq,i^)(bl,,rKu,T+ \)- (19) 



q,VT 



Here io-j-ik, u) is the frequency of the normal lattice vibrations of the mode r (r = 1, 2, . . . , 12) 
with the longitudinal wavenumber k and the azimuthal quantum number v. The adimensional 
wavenumber (quasi- momentum) along the nanotube, k = x^i' t^'^es quasi-continuum values 
(for L S> 1) in the range — tt < k < tt. The azimuthal quantum number takes discrete values 
u = with n2 = 0, ±1, . . . , ±^ if N is odd and n2 = 0, ±1, . . . , ±(f - 1), f if is even. 

The frequencies LOrikjv) and the coefficients of the transformation (17) can be found from 
the diagonalization condition of the potential energy of the lattice displacements in (16) with 
the orthonormalization conditions 

Y2 E (Krik, l^)U,,r'{k, U) + Sl,{k, u)S,,r'{K v) + VlAK ^)V,,Ak, ^)) = Kr', 
Q 

1 UlAk, iy)U^Ak, ^) = J2^ Sl^ik, u)S^,r{k. ^) = ]^ E Krik. ^)^.',r(/c, 1^) = 5,,^, 

T T T 

E^.V(^'^)^^?V(^'^) =E^^,r(^,^)V;',r(fc,i^) =E^.:r(^'^)^.',r(^,^) =0. (20) 
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Note that any linear form of lattice displacements, such as WSm,n,g, ^Sm,n,e and Cm,n,Q-, 
after applying the transformation (17) can be written as 

^"^'"■^ = ~7mWz ]^ ^\'')Q^,'^,r (21) 

where FQ{k,i'\T) is a linear form of the transformation coefficients Sg^r{k,i^), V^,T(fc,z^) and 

Therefore, in general, the frequencies a;,- (A;, v) of the normal vibrations of the lattice can be 
represented as 

4{k, ^) = 1 ^( A J2[\W6,{k, Hr)|2 + \m,{k, Hr)|2] + ^|C,(fc, i^\t)\^) . (22) 
e s 

The equations for the normal modes of lattice vibrations are too complicated to be solved 
analytically in the general case. For carbon nanotiibcs the phonon modes were calculated numer- 
ically (see, e.g., [1, 28] and references therein). Here we do not calculate the phonon spectrum 
explicitely, instead we use the general relations (20), (22) to get some estimates which only de- 
pend on the parameters k, kc and M. Meanwhile, the explicit expressions for the electron 
dispersions are more important for us and will be derived below. 



4 Electron Hamiltonian 

The electron eigenstates are found from the tight-binding model using the nearest-neighbour 
hopping approximation. In this approximation the Hamiltonian which describes electron states 
is given by 

ee,a S 

Here a|g^ (033,0-) are creation (annihilation) operators of a 7r-electron with the spin a on the site 
86, £0 is the vr-electron energy, J is the energy of the hopping interaction between the nearest 
neighbours and the summation over 6 denotes the summation over the three nearest neighbour 
sites. 

By the unitary transformation 

am,n,g,a = E e'''"+^''^u,,,{k, v)ck,u,X,., (24) 

^ k,v,X 

with ^ 

4 E^e,A(^''')v(^''^) = <^A,A' (25) 
Q 

the Hamiltonian (23) is transformed into a diagonal form (see Appendix 1): 

He= ^>^ik^'')4,u,X,a(^k,u,X,a- (26) 

k,i/,X,a 

Here k is an adimensional quasi-momentum along the nanotube, v is an azimuthal quantum 
number, and A = 1,2,3,4 labels the four series (due to the four atoms in each cell), of ID 
electronic bands with the dispersion laws 

Ex{k,i^) = So ±£±ik,iy), (27) 
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where 



£±ik,u) = JWi + 4cos2(^)±4cos(^)cos(^) 



(28) 



In (26) the operators c|, (^{ck,u,\,<7 ) are creation (annihilation) operators of electrons with the 
quasi-momentum k and spin a in the band {v, A). If we label the electronic bands as 



Ei{k,u) = So - £+{k,u), E2{k,u) = Eq - S-{k,u), 
E3{k,iy) = So + S-{k,u), Ei{k,v) = Eq + S+{k,u), 

then the matrix of the unitary transformation coefficients u (24) is given by 



(29) 



u{k, u) = 



r k — iy 



_^-i(!^-e+) 



ii!^-e+) _^i{h^+e.) 



(30) 



where the phases satisfy the relation (148), given in Appendix 1. 

5 Electron-Phonon Hamiltonian 

The electron-phonon interaction originates from different mechanisms [10, 14, 11, 7, 8]. Usually, 

the dependence of the hopping interaction between the nearest neighbours J(ac);5(ac) o^i their 
separation is considered and in the linear approximation with respect to the displacements one 
has 



J(.. 



J - G2W5^ 



(31) 



'(s);<5(s) 

In general, neighbouring atoms also alter the energy of the 7r-electrons on each site and so, in 
the same linear approximation, we can write 
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(32) 



Thus, the total electron-phonon interaction Hamiltonian should be taken in the following 
form 

Hint = ( «L,(7"ae,(T [Xl X] + ^2 C'as] + ^2 «le,cT«5( 

W5^), (33) 

where we have used the translation index operator 5(«) defined in (8). 

The unitary transformations (24) and (17), transform the interaction Hamiltonian into 



Hint — 



2VSLN 



k,u,X,X' ,q,^,T,cr 

where Qq,fi,T was determined in (18) and 
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8^,8 



TeAl^ Mk) = Wg{q, ii\t) + C^(g, //|t). 



(34) 



(35) 



Note that Tg>^Q{k,i^;q, h\t) = T*^^,{k + q,i^ + n; -q, -fi\T) and that Ti^3 = T^^i = T2,a 
^4,2 = 0. The diagonal elements, at q' = g, are 



(36) 
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and the nonzero ofF-diagonal elements, Q q' ■, are given by 



G2 

Te',g{k, v; q, /x|r) = -y= Wg>^g{k, u; q, /x|r), (37) 
s/ M 



where Wg{q, /^Ir), Cg{q, /x|r) and Wg'^g{k, u; q, ij,\t) are determined only by the coefficients of the 
phonon unitary transformation (17). In particular, 

Wiiq,,,\r) = V3sin(|)(^Si,, + e-^tcos(|)52,0) + 

+ iV3cos(|) sin(|)e-'t^2,r + cos(|)e-*tF2,, - e-'^l^4,r, 

W2{q,^i\r) = V3sin(|)(52,. + cos(|)e^t5i,0) + 

+ zV3cos(^) sin(|)e^^?7i,^ - cos(|)e^^yi,^ + Fs.r, 

W3{q,ix\T) = V3sin(^)(53,. + cos(|)e^t54,.)) + 

+ zV3cos(^) sin(|)e^tc/4^^ + cos(|)e^ty4^^ - V2,r, 

W4{q,l^\T) = ^/3sin(|) (^54,. + cos(|)e-*t53,^)) + 

+ zV3cos(|) sin(|)e-tc/3,, - cos{^)e-'^Vs,r + e'^Fi,,, (38) 



Wi2{iy;q,fi\T) = e-t(V3sin(^)(^cos(^)Si,, + e-^tcos(^)S2,.))- 

V3cos(^) (smQu,,r - e-'t sin(^)t/2,,; 



and 



+ e-tcos(^)F2,.-cos(^)Fi,.), 
T^i4(A;;g,//|r) = e'"' (Vi,r - e-'iV^^r) , W23{q, ii\t) = V3,r - V2,r, 
Ws4{'^;q,fi\T) = e*f(V3sin(|)(cos(^)53,, + e^tcos(^)54,.)) - 

- iv^cos(^) (sin(^)C/3,, - e^^ sin(^)C/4,,)) + 

+ e^tcos(^)T4,r-cos(^)F3,r), (39) 

where Sg,r = Sg,r{q,IJ-), Vg^r = Vg^riQ.lj), and Ug^r = Ug^Aq,!^-)- 

Thus, the functions F^y{k,q;u, /j,) are determined by the interaction parameters xi, X2, 
G2 and the coefficients of the unitary transformations (24) and (17). Note that, in (34), the 
azimuthal numbers satisfy the relation vi = i^+ii, for which the following rule should be applied: 
if + > 27r, then ^ v[ = vi ± 27r in such a way that < 27r. 

6 Adiabatic approximation 

The total Hamiltonian of the system is then given by 

H = He + Hph + Hint, (40) 
where H^, Hp^ and Hint are given by (26), (19) and (34), respectively. 
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Below we consider only one-particle states in a carbon nanotube taking into account the 
interaction of the particle with the lattice distortions. When the couphng constant of this 
interaction is strong enough, this interaction can lead to the self-trapping of the particle. The 
self-trapped states are usually described in the adiabatic approximation. In this approximation 
the wavefunction of the system is represented as 

|^)=C/|Ve), (41) 

where U is a unitary operator of the coherent atom displacements induced by the presence of 
the quasiparticle and so is given by an expression of the form 



U = exp 



(42) 



and | (/'e) is the wavefunction of the quasiparticle itself. Moreover we require that it satisfies 

(V'elV'e) = 1. 

In (42) the coefficients Pf^^riQ) depend on the state of the quasiparticle which, in turn, is 
determined by the lattice configuration. Using (41) in the Schrodinger equation H\^) = 
we find the equation for the electronic part Itpe) of the total function (41) 

H\^e) = E\4^e), (43) 

where 

H = U^HU = W + He + Hint + Hph + Hd. (44) 

Here 

W nur{q,lj^)\l3r{q,l^t (45) 

is the energy of the lattice deformation, 

k,V,\,(T 

k,q,\,X ,T,(7 

is the Hamiltonian of quasiparticles in the deformed lattice with the deformation potential given 
by 

Qriq, /^) = (2^;;;^) ' (Mq, m) + P*A-q, -f^)) , (47) 

and 

Hd= ^^T{q,fi){Mq,l^)bU,r + P*riq,f^)bq,i^,T) (48) 

is the displacement term of the phonon Hamiltonian. The latter term, H^, is linear with respect 
to the phonon operators and appears here as a result of the action of the unitary operator (42). 
With the help of the unitary transformation 

Ck,u,\,a = Yi'q;\{k,y)Cr,,a, (49) 

we can introduce the new Fermi operators C^^a for which, in the general case, the quantum 
number r/ is a multicomponent index. The coefficients tpn-xik^u) are to be chosen from the 
condition that the electron Hamiltonian (46) can be transformed into a diagonal form: 

He = J2E,Cl,Cr,,a. (50) 
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This requirement leads to the following equations for the transformation coefficients: 



+ o /^TTT E\\'{k-q,q\v- IJ.,iJ)QT{(l,IJ')'il^r,-,\'{k-q,v- ii). (51) 

Solutions of this system of equations, with the orthonormalization condition 

^)'4^v'-Ak, v) = Sr,,r,' (52) 

then give us the coefficients ipri;\{k, v) as well as the eigenvalues of the electron energy levels. 
After the transformation (49) the interaction Hamiltonian becomes 



E ^li'i<l^t^)Cl.C,,^„Q,,,,r, (53) 



where 



= E rr,Mk + Q^^ + f^)Px'!xik,<l-,'^,t^)^v'-Ak,^)- (54) 

The operator Hint can be separated into two parts. The most important term, Had, is 
the diagonal part of Hint with respect to the electron quantum numbers r] [r] = r]' in (53)). 
The remainder, Hna, the off-diagonal part of Hint, corresponds to phonon induced transitions 
between the adiabatic terms determined by Eqs. (51). So we can represent the Hamiltonian 
(44) in the form 

H = Ho + Hna (55) 

where 

Ho = W + He + Had + Hph + Ha (56) 

describes the system in the adiabatic approximation and Hna is the nonadiabaticity operator. 
At large enough electron-phonon coupling, the nonadiabaticity is less important and the 

operator Hna can be considered as a perturbation. In the zero-order adiabatic approximation 
the quasiparticle wavcfunction \ipe^^) does not depend on phonon variables. In the case of a 
system with Ng electrons it can be represented as a product of A^e electron creation operators 
which act on the quasiparticle vacuum state. 

In particular, the one-particle states are described by the function 

lV'f)=Ct,JO), (57) 

where |0) is the quasiparticle and phonons' vacuum state, the index r] labels the adiabatic 
state which is occupied by the quasiparticle. For the ground state we put rj = g. The total 
wavefunction of the system (41) describes the self-trapped states of a large polaron in the zero- 
order adiabatic approximation. 

Note that the function (57) is an eigenstate of the zero-order adiabatic Hamiltonian Hq: 



Ho\^f^) = \W + Eg 



+ {ihu;Aq,fi)P,A<l) + W 6LArat(g,/x) ^^^*^^'^^^^^-^-- + ^•'•)] '^^"^^ ^^^^ 
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with the energy Sg = W + Eg provided that the coefficients /3r (<?) in (42) satisfy: 



The adiabatic electron states are determined by (51) in which the lattice distortion Qt{c[,IJ'), 
according to Eqs. (47,59), is self-consistently determined by the electron state: 

= -^Vm 'a;?(g,/x) '^^^^^^^^ + ^' " + 

Substituting (60) into equations (51) for the occupied electron state, we obtain a nonlinear 
equation for xjjg-x^k,^) whose solution, satisfying the normalization condition (52), gives the 
wavefunction and eigenenergy Eg of the electron ground state and, therefore, the self-consistent 
lattice distortion. All other unoccupied excited electron states with r] ^ g can be found from 
the linear equations (51) with the given deformational potential. 

Using the inverse unitary transformations (49) and (24), we can rewrite the eigenfunction 
(57) in the following form: 

= E V'.;A(fc,^)i,,,,jo) = ^ v.,.4,Jo), (61) 

\,v,k ae 

where 

= ITTfTr E e^^^™+'^")«,,A(fc, v). (62) 



Here "^g,^ is the polaron wave function, i.e., the probability amplitude of the distribution of a 



quasiparticle over the nanotube sites: P(ae) = |V'g 



7 Large polaron state 

Putting Eqs (60) into (51) gives us the nonlinear equations 

(Ey,{k,v) - E)i,y,{k,v) = 

E Gi^:i"{ \ l\[ Ov'^,(fci,^i)VAi(fci+9,^i+/x)V'A'(fc-9,^-/^) (63) 

A',Aj,Ai,fei,i/i,q,/i V ' ' / 

for the one-electron ground state. Here, and from now onwards, we have omitted the index 
rj = g and introduced the notation 

Al,Ai ( k, ku q \ _ I ^ ^xM'' -q,i^-fi;q, fi)F^^]l{kuiyi;q, ^i) 
^A,A' M j"l2^ u:l{q,n) ■ ^^^^ 

We see that all sub-levels of all sub-bands participate in the formation of the self-trapped 
electron states and, in general, there are many solutions of Eq.(63). Among these solutions there 
are 'one-band' solutions in which only the function V'a(^) ^) with quantum numbers A = Aq and 
v = vq IS nonzero and all other functions ip\{k,i') = with A 7^ Aq and v ^ vq. But not all of 
these solutions are stable. 
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Next we consider the 'one-band' self-trapped state which is stable and is split off from the 
lowest energy subband in (29), namely from £'i(A;, 0) with Aq = 1 and uq = 0. In this case Eq. 
(63) becomes 



= - E^{k,(}))il;i{k,0) 
1 



+ 7^E<^(^'^i''?)V'Kfei>0)Vi(^i+9,0)V'i(fc-g,0), (65) 

ki,q 



where 



G(fe,fei,g) = Gy ( J' ly (66) 



To solve (65), we introduce the function 

¥'(O = ^Ee''"^i(^'0) (67) 

which depends on the continuous variable that is a dimensionless coordinate along the nan- 
otube axis related to z hj ( = z/3d. 

Then we assume that in the site representation a solution of (65) is given by a wave packet 
broad enough so that it is sufficiently narrow in the k- representation. This means that ijji{k, 0) 
is essentially nonzero only in a small region of fe-values in the vicinity of = 0. Therefore, we 
can use the long-wave approximation 



Ei{k,0) = So-J^5 + 4cos{'^)^Ei{0) + ^Jk^, 
G{k,ki,q) « Go(0,0,0)=G, (68) 

where 

Ei{0) = So - 3J (69) 

is the energy bottom of the subband Ei(k,0). 

Using Eqs. (35) - (39) and (19) for = = in the long-wave approximation, we can 
represent the nonlinearity parameter G as 

where ai is a constant of the order of unity, while the constants 61,62,0 are less than 1. Intro- 
ducing 

A=E - Ei{0), (71) 
we can transform Eq.(65) into a differential equation for ip{C). 



A^(C) + ^^ + |l^(C)lV(C) = 0, (72) 



which is the well-known stationary nonlinear Schrodinger equation (NLSE). Its normalized so- 
lution is given by 

= Vfcosh(,o(c-Co)) ^''^ 

with the eigenvalue 

Ao = (74) 
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where 

.0 = §. (T5) 

Thus, the eigenenergy of this state is 

^° = ^^^°)-4JiV^- ^^^^ 

The probabiUty amphtude (62) of a quasiparticle distribution over the nanotube sites, in 
this state, is given by 

^m,„,, = — ^Ee*'"^^,i(^,0)Vi(fc,0). (77) 



2VLN , 

The cxpUcit expressions for Ug_i{k,0) are given in (30). In the long- wave approximation for the 
phase 9-^-{k, 0) we find from (148) that 0^{k, 0) k/l2. Then, using the expressions for Ug^i{k, 0) 
and taking into account the definition (67), we obtain 

i'm,n,e = ■T^fizm,e), (78) 



where Zm,Q are the atom positions along the nanotube axis (2): 

Zm,i = {m- ^)3d, Zm,2 = (jn - ^)3d, 

Zm,3 = (m + ^)3d, Zm,4. = {m + ^)3d. (79) 

Therefore, according to our solution (73), the probability distribution of a quasiparticle over 
the nanotube sites is given by 

Thus, the quasiparticle is localised along the tube axis and distributed uniformly over the tube 
azimuthal angle. Therefore, (80) describes a quasi-lD large polaron. In this state, as well as in 
other one-band states, according to (60), only the total symmetrical distortion of the nanotube 
takes place, i.e. Qr('?iO) / with jj, = and QriQ^fJ-) = for /Lt 7^ 0. The total energy of the 
polaron state according to (58), is 

So = W + Eo = E,{0)-^, (81) 
and, thus, depends on the diameter of the nanotube. 



8 Transition to states with broken Eixial symmetry 

As we see from (78), (73) and (80), our solution, obtained in the long- wave (continuum) approx- 
imation, possesses the azimuthal symmetry and describes a quasi-lD large polaron state which 
is localized along the nanotube axis in the region Az = Moreover, (75) shows that as the 
electron-phonon coupling increases the region of the localization gets smaller. Consequently, the 
wave packet in the quasimomentum representation becomes broader and the electron states with 
higher energies participate in the formation of the polaron state. At strong enough coupling the 
long-wave (continuum) approximation is not valid. Moreover, the electron states from the upper 
bands can also contribute to the polaron formation. To consider the transition from the large 
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polaron state to the small one, it is convenient to transform Eqs.(63) into the site representation. 
As a first step, let us introduce the functions 

4>e{k, J^) = 7^J2 ^e,A(^) v). (82) 

Then Eqs.(63) can be rewritten as the following system of equations 
E(f)i{k,u) = SQ(t)i{k,v)-2Jcos{^)e-''i^2{k,u)- Je-^^^i{k,v)- 

ki,i'i,q,ij.e'i,ei,e' 

E(j)2{k,iy) = So(l)2ik,u) -2Jcos(^)e^'^(t)i{k,v) - J(l)s{k,v) - 

~ Yn ^ ^ T^^f\k,iy]ki,i^i-q,ii)(pl^iki,i^i)(pg'^{ki+q,iyi+iJ,)(pg>{k-q,i^- IJ,), 

E(t)^{k,v) = £Q(l)-i{k,v) -2J cos{^)e'-^(t)i{k,v) - J(l)2{k,v) - 

V _-v i. 

E(f)i{k,iy) = £o(t)4{k,iy) - 2J cos{-)e '2 03(/c, z^) - Je z^) - 

~ Yn Yl Yl '^4^fHk,J^\ki,ui;q,ii)(t)*^{ki,ui)(f>^,^{ki+q,ui+ii)(f>g,{k-q,u- 



(83) 



where 



r^;^, {k,u;k,,ur,q,n) = -^ . (84) 

In the derivation of these equations we have used the explicit expressions (64) and (35), the 
orthonormalization conditions (52) and the following expressions for £±{k,i') 

£±{k,u) = J (^2cos(^) ±e-*i) e±2*e±(fc,v) ^ j (^2cos(^) ± e^t^ e^2ie±(fc,v)_ ^g^^ 



To describe the system in the site representation, we introduce 

1 
L 



^g,m{^) = jY^'''^Mk,'^), (86) 
k 



and obtain 

2 

E(t)l,m{l^) = So(f)l,m{l^) - 2J COs{^)e-'^(t)2,mM - J^4,m-l(z^) " ^7 X) {^Xl^l,m{^l)^l,m{j^l + + 
+ XlG2[cOs(y)e*^^2,m('^l)?!'l,m(j^l + /^) + 

+ C0S( —)e '~0t,„(z^l)02,m(i^l +1^)+ (t)*l,m{^l)(t^4,m-l{l^l + 1^) + 

+ (Pl,m-l{^l)<Pl,7n{l^l + l^)])4'l,m{k - q,V - ix) - 
O? ( 

~ nY [XlG2[(t>*l,m{'^l)(t>l,mi'^l +fJ-) + 
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+ cos( ^^ ^ ^ )e~'^(/>i^^(z/i)02,m(t^i + ^)])e"*^(/)2,m(fc -q,v- ^i)- 
o? / 

+ 2G2[<?!)4_m_i(l^l)(/'l,m(j^l +//) + 4>l^^{vi)(l)i^ra-l{vi + lj)'\)<t>i,m-l{k -q,l^~ jj)- 

with similar equations for ^> = 2, 3, 4. 

When deriving equations (87) we have made a qualitative estimate of the expressions of the 
form 

12^^ ^^^(^ ^''^ 

by assuming that the main contribution to these quantities comes from the lattice variables with 
small q and fi. This gives us an estimate of a\ in Eq.(87). 

In zigzag nanotubes, one can identify zigzag chains of carbon atoms which encircle the 
nanotube. Let the ring chain j consists of atoms enumerated as (m, n, 1) and (m, n,2), the 
zigzag chain of atoms (m, n, 3) and (m,n, 4) corresponds to the ring j + 1, and the chain of 
atoms (m — 1, n, 3) and (m — 1, n, 4), respectively, to the ring with number j — 1. Then we can 
enumerate atoms as {j,n,p) where p = 0,1. Note that the indices {j,p) coincide with the ones 
used in the numerical calculations [23]. 

A circle around the nanotube is a zigzag ring chain, with two atoms per unit cell, which 
contains 2N atoms. The atoms of the j-th chain are equivalent except that atoms with p = 
are coupled to the (j — l)-th chain and those with p = 1 to the {j + l)-th chain, and these two 
sets of atoms are shifted with respect to each other in the opposite directions from some central 
line, Zj, (symmetry axes). Thus, we can put 

1 1 /^"^ -^^^ \ 



. 2N-1 

^ E (89) 



\n=l n=l 



V2N 

We see that il^jiv) is a /^-representation for a simple chain with 2N atoms with k = u/2 = 

imi/n. Therefore, this zigzag ring chain is equivalent to an isolated nanocirclc, studied in [25]. 

Introducing the notation: 4>i^m{y) = 4^o,j{y)i 4>2,m{^) = fpiji^) and neglecting X2 we 

can rewrite Eq.(87) as follows: 

E(l>o,j{i^) = - 2J cos(-)<^ij(i^) - J0ij_i(i/) - 1^ E '/'o,j(^i)</'o,j(i^i + ^J)<l>o,j{^ - f^), 

U1,IJ, 

(90) 

where G is given by Eq. (70). 

For the azimuthal symmetric solution the only nonzero functions are those with zero argu- 
ment, 1^ = 0. In this case we can use the continuum approximation: 

<f>o,j = 0(Co,j), <f>i,j = 0(Co,j + \) = HCoj) + ^0'(Coj) + ^<l>"{Co,j), 

<pi,j-i = <p{Co,j-l) = m,j)-l^'iCo,j) + ^(^"{Co,j)- (91) 

As a result Eq. (90) transforms into the continuum NLSE (72). The azimuthally symmetric 
solution of this equation does not always correspond to the state of the lowest energy. To find 



16 



the lowest energy state, we consider Eq. (90) assuming that the electron is localized mainly on 
one chain (for simplicity we label it by j = 0) and we look for a solution of the form 

cPpjiiy)=Ap^,^{iy), (92) 

where Apj are given by Eq. (73) with (^q describing the position of the considered chain. 

We can now consider the equation for the chain j = only. For 4){i') we obtain the equation: 

{E - £{u))4>{u) + _^^0*(i.i)</,(i.i+^)</>(i.-^), (93) 

where 

£{u)=£o- -iJcos"^. (94) 
^0,0 2 

Moreover, from (73) we find that 

A A [so 1 A [go 1 (nr.\ 

^0,0 = ^1,0 = \ H:, TT 7T7^) M-l — \ — rrr 77^. ydO) 

V 2 cosh(5(o/12) V 2 cosh(55o/12) 

Assuming that the function </>(z^) is essentially nonzero in the vicinity of the zero values of 
V, the energy dispersion can be written in the long-wave approximation as 



where 



£{k) = £{0) + J( 2 ) +•••> (96) 

£{0) = £,-j(2 + ^ j . (97) 
To solve Eq.(93) we introduce the function 

</^(^) = ^E«'''^W (98) 

of the continuum variable x with k = v/2 being the quasimomentum of the nanotube circle with 
one atom per unit cell. Therefore, the quasimomentum representation should be studied in the 
extended band scheme and — tt < k < tt. Note that f{x) is a periodic function, ip(x+2N) = ^p{x), 
and that the discrete values of x = n, n = 1,2, ...2N — 1, correspond to the atom positions in 
the zigzag ring. 

Using the approximation (96) one can transform (93) into a nonlinear differential equation 
for ip{x) (stationary NLSE): 

+ GAl,\<p{x)\M^) + A^ix) = 0, (99) 

where A = £;-f(0). 

As it has been shown in [25], the solution of Eq. (93), satisfying the normalization condition 

/ \cp{x)fdx = l, (100) 
Jo 

is expressed via the elliptic Jacobi functions: 

'2K{k)x 



2N 



-,k 



(101) 
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Here g = GA^ q/ {2J), and K(fc) and E(fc) are complete elliptic integrals of the first and second 
kind, respectively [29]. The modulus of the elliptic Jacobi function, fc, is determined from the 
relation 

E(fc)K(fc) = (102) 



2 



The eigenvalue of the solution (101) is 



According to [25], the azimuthally symmetric solution exists (relation (102) admits solution) 
only when g exceeds the critical value of the nonlinearity constant 

9> 9cr = (104) 



or, in an explicit form, 

At ^ > 1, (105) 

27r2 cosh2(a/(4A^)) ' ^ ' 

where a = G/J is the adimensional electron-phonon coupling constant. 

Prom (105) we find the critical value of the coupling constant, acr = 2.6 for N = 8. According 
to the numerical solution [23], the critical value of the coupling constant (xi + G2)^ / {kJ) ~ 3.2 
for this value of N . Comparing this result with the analytical prediction, we conclude, that the 
parameter oi in (70) is ai k, 0.9. Therefore, the estimation ai 1 made above is indeed valid, 
which justifies our analytical results. Of course, the applicability of this approach far from the 
transition breaks down because the continuum approximation itself is not valid anymore. 



9 Large polaron states in semiconducting nanotubes 

In zigzag nanotubes, when TV is not a multiple of 3, there is an energy gap in the electron 
spectrum (29). In the carbon SWNT this energy gap, A, separates the ID electron sub-bands 
£q — £±{k, v), which are fully occupied, from the empty ones with energy £q + £±{k, v). Such 
nanotubes are semiconducting [1]. Their charge carriers are either electrons (in the conducting 
band) or holes (in the valance band) . For semiconducting zigzag nanotubes N can be represented 
as N = 3no + 1 ov N = 3no — 1. The lowest conducting subband above the energy gap is 

E3{k,uo) = £o + S^{k,uo), (106) 
and the highest valence subband below the gap is 

E2{k,uo) = So - S-{k,uo) (107) 
with uq = 27rno/-/V. So, the energy gap in semiconducting carbon nanotubes is given by 

A = E3{0,uo) - E2iO,iyo) = 2f-(0,z.o) = 2J|1 - 2cos(^)|. (108) 

Next, we consider a self-trapped state of electrons in the lowest conducting band E2,{k,v). 
Because £±{k,—i') = £±{k,i'), this subband is doubly degenerate. In this case we look for a 
solution in which only the functions ipsi^, ^) with v = ±z/o are nonzero. Then Eqs. (63) become 

E\l)3{k,v) = E2{k,uo)ilj3{k,u) - 

k',q, v' 

+ Gt^l,(fe, k\ q) rsik', -iy)Mk' + q, ^)Mk - Q, -i^)) , (109) 
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where v, v' = ztz/Q and 

= J/ 2). (110) 

Gi,!U*,fc'..) = (111) 

Here the equivalence of the azimuthal numbers ^ and it 27r should be taken into account. 
To solve (109), we introduce functions of the continuum variable x using the relation (67) 

<^.,3(x) = 4f Ee'''^3(fc,^) (112) 

and use the long-wave approximation 

4y(^>^''9) ~ G(;i,(0,0,0) = Gi, 

G^J^_,{k,k',q) « g(^L(0,0,0) = G2. (113) 

Note that 

E3iO,uo) = £0 + (114) 
is the energy bottom of the subband £^3(fc, fo) and 

211 -2cos(^)|^2 , , 

is the quasiparticle effective mass in the subband E^lkjUo). 

In this case Eqs.(109) are transformed into a differential equations for (p,^^3{x): 



A¥'.,3(x) + ^^-^0^ + ^(Gi|¥'.,3(x)|2 + (Gi + G2)|^-.,3(x)|2)^.,3(x) = 0, (116) 



where 

A =E - ^3(0,1/0), (117) 

and u = it'o- 

We observe that equations (116) admit two types of soliton-like ground state solutions. The 
first type corresponds to the given azimuthal quantum number: v = vq 01 v = —vq. Such 
solutions describe solitons with a fixed value of the azimuthal number and are formed by the 
electron sublevels with 1/ from the doubly degenerate band, i.e., only one function ip^, ^ is 
nonzero and the other one vanishes: ip-i, = 0. These types of solitons are described by the 
NLSE: 

t2 J2.„ ^(^\ ^ 



A¥'.,3(x) + — ' + ^|^.,3(a:)|V.,3(^) = 0. (118) 
A normalised solution of this equation is given by 

^Z,u{x) = T (119) 

V 2 cosh(5ia;) 

with the eigenvalue 

A.^-|A, (1.0) 
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where 



_ mGi _ Gi|l-2cos(^)| 
^'~2h'N~ JiV cos(^) • 



Thus, the eigenenergy of these states is 

1 G?|l-2cos(^)| ^ ^ 

^^=^°+2^- 4JiV^cos(t) • ^'^'^ 

The energy levels of the two solitons with different azimuthal numbers {v = vq and v = —vq) 
arc degenerate, similarly to the levels of the corresponding bands. However, according to Jan- 
Teller theorem, this degeneracy can be broken by the distortions of the lattice resulting in the 
hybridization of these two states. 

Next we consider a case when both these functions are nonzero, ^± 7^ 0. In this case ip± are 
determined by the system of nonlinear equations (116). A normalised solution of these equations 
is given by 

^±vo = ■^e''^" 95/1,3, (123) 

where (f)± arc arbitrary phases and where 99/1 3 satisfies the NLSE (118) in which the nonlinearity 
parameter Gi is replaced by G\ — > (2Gi + G2)/2. Therefore, this solution is given by (119) 
with 

m{ 2Gi + G2) _ (2Gi+G2)|l-2c os(^) 
A-h^N ~ 2JNcosif] 



Its eigenenergy is 



p , 1. (2Gi + G2)^|l-2cos(^)| 
E,-So + -A i6Jiv2eos(f) " ^^^^^ 

This hybrid soliton possesses a zero azimuthal number while solitons (119) have a nonvan- 

ishing one: i' = uq or u = —uq. The energy level of the hybrid soliton state. Eh, is lower than 
the level of a soliton with the fixed azimuthal number, Ei, because {2Gi + G2)/2 > Gi. Note 
also that the deformation of the nanotube in this state is more complicated due to the fact 
that the components Q±2uo of the lattice distortion, as well as the Qo-component, are non-zero. 
Moreover, the probability distributions of a quasiparticle over the nanotube sites in the state of 
a hybrid polaron and in the state of polaron with a fixed azimuthal number, are different. 

For a polaron state with a fixed azimuthal quantum number (e.g. f = vq), the probability 
amplitude (62) is 

V^^,„,, = — ^ e^^'"^+'^°")n,,3(A;, yo). (126) 



2^LN , 

The explicit expressions for uj^^^k, v) are given in (30). In the long-wave approximation for the 
phase 6-{k,i'o) we find from (148) that 

'"^''-'*'°'" 4(2cob('h.)-1) - 
Then, using the expressions for Uj^3{k,i') and taking into account the definition (67), we obtain 

^m,n,i = ^e-°("-i)(^,,3(m-i + ^<5), 

V'm,n,2 = --^e-o(-+i)^.,3(m-l-i<5). 



^m,n,3 = --^e-o('^+l)v..,3(m+^ + ^<5). 



V'm,n,4 = -^e-°("-i)<^.,3(m+i-i5), (128) 
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where 

Therefore, according to (128) and (119), the quasiparticle distribution over the nanotube 
sites in this state is 

PM.n) = = ^-jj^^^^-i-^^, (130) 

where Zm,e are the atom positions along the nanotube axis (2), the "+" and "— " signs correspond 
respectively to atoms with an odd {g = 1,3) and even {g = 2,4) index g . Usually these two 
types of atoms are labelled as A and B atoms. 

We see that the quasiparticle is localised along the tube axis and is uniformly distributed over 
the tube azimuthal angle like a quasi-lD large polaron. But the distributions of the quasiparticle 
among A and B sites are shifted relatively to each other by the value 3d5{uQ). 

For a hybrid polaron state which possesses zero azimuthal number, the probability amplitudes 
(62) are 

cos (i^o(n - \) + (t)oj 1 1 

^m,n,l = -7^ </^m(^- 3 + 2^)' 



cos(z/o(n+i) + ^o 



1 1 



cos (^vo{n+ i) + 00 



1 1 



^m,n,3 = ¥'m("^+ g + 2*^)' 

COS (i^o(n - I) + (t)o] 1 1 

^m,n,4 = ^^1^ </'m(^+ 3 - 2^)' 

Therefore, the quasiparticle distribution over the nanotube sites in this state is given by 

P(mn)-P(z d> )- 9hCosHnf,cPn,e + 4>(i) ^32^ 

4Ar cosh [^{zm^g ± i^do)) 

where (pn.Q is the angle for the atoms position in the nanotube (2), (jj^^g = na for ^ = 1,4 
and (pn,^ = (n + for ^ = 2,3; no is a number which determines the azimuthal number vq 
{i/Q = 2Trno/N), and the " + " and " — " signs correspond to the odd and even values of £> as 
above. 

We see, that in this polaron state the quasiparticle is localised along the tube axis and is 
modulated over the tube azimuthal angle with the angle modulation 27r/no. The longitudinal 
distributions of the quasiparticle among A and B sites are shifted relatively to each other by 
the value 3dS{vo). 



10 Conclusions 

In this paper we have derived the equations describing self-trapped states in zigzag nanotubes 
taking into account the electron-phonon coupling. Wc defined the electron and phonon Hamilto- 
nians using the tight-binding model and derived the clcctron-phonon interaction arising due to 
the dependence of both the on-site and hopping interaction energies on lattice distortions. Next 
we performed the adiabatic approximation and we obtained the basic equations of the model. 
These are the equations in the site rcpresantation that were used by us in [23] to compute nu- 
merical solutions of nanotubes states and to determine the ranges of parameters for which the 
lowest states were soliton or polaron in nature. 
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In this paper we have studied this problem analytically. We have shown that the electrons in 
low lying states of the electronic Hamiltonian form polaron-like states. We have also looked at 
the sets of parameters for which the continuum approximation holds and the system is described 
by the nonlinear Schrodinger equation. This has given us good approximate solutions of the full 
equations (thus giving us good starting configurations for numerical simulations) and has also 
allowed us to compare our predictions with the numerical results [23]. 

Our results demonstrate the richness of the spectrum of polaron states. They include quasi- 
ID states with azimuthal symmetry for not too strong coupling constant, and, at relatively high 
coupling, states with broken azimuthal symmetry which are spread in more than one dimension. 
Theoretical estimates of the critical value of the coupling constants between the two regimes 
of self-trapping (with or without axial symmetry) are in good agreement with our numerical 
results [23]. 

We have also found that for the values of the parameters corresponding to carbon nanotubes, 
the lowest energy states are ring-like in nature with their profiles resembling a NLS soliton, ie 
similar to a Davydov soliton as was claimed in [23]. 

We have considered the polaron state of an electron (or a hole) in semiconducting carbon 
nanotubes and have shown that the degeneracy of the conducting (or valence) band with respect 
to the azimuthal quantum number plays an important role. The polarons with lowest energy 
spontaneously break down the azimuthal symmetry as well as the translational one and posses 
an inner structure: they are self-trapped along the nanotube axis and are modulated around the 
nanotube. 

Next wc plan to look in more detail at some higher lying states and study their properties. 
We are also planning to study the electric conduction properties of our system. 
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12 Appendix 1. Diagonalization of the polaron Hamiltonian 

Due to the fact that the diagonal expression 

He,0 = £oYl «t<Taae,a (133) 
ae,c7 

remains diagonal under any unitary transformation, wc consider only Hj. Omitting the multi- 
plier J and the spin index we note that Hj is given by, in an explicit form , 

Hj = — l^g]^ „j(am,n-1.2 + Qm.n.2 + am-hnA) 

m,n 

+<^iL,n,2('^m,n,l + am,n,3 + am,n+l,l) 
+"m,n,3('^»n,">4 + «m,n+l,4 + am,n,2) 

+"m,n,4('^m,n-l,3 + CLm+l,n,l + Om.n.s) • (134) 

Due to the translational invariance (with respect to shifting the index m) and the rotational 
invariance (changing n) we can perform the transformation 



' KM 



=i:e^'"'^''-c^k,u,„ (135) 
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which transforms the Hamiltonian (134) into a diagonal form with respect to the indices k and 
v and we get 



+ 2 cos ^ (e^t al^^^ak,u,4 + e~^^ 4 ,, 40^,^,3) + 



+ 2cos|^(e '^'4,.iafc,i.,2 + e*2al 2«fe,i',i 



Note that a direct way to diagonahse (136) is to use the unitary transformation 



with 



(136) 

(137) 
(138) 



A j 
which leads to a system of equations (four in our case) for the coefHcients UQ^\{k^u) which 
diagonalisc the Hamilionian. A solution of these equations would give us the coefficients of the 
transformation as well as the eigenvalues Ex{k,u) (A = 1,2,3,4). Instead, we prefer to use a 
sequential diagonalization. 

To do this we choose any two different pairs of operators 0^,1/,^ and a^^j^^gi and using unitary 
transformations we first diagonalise two of the four lines in (136). Taking the following two 
pairs: {ak,u,i, 0'k,u,2} and {ak,u,3, o,k,i',4}, and diagonalising the last two lines in (136) by the 
unitary transformations, we get 



and 



0-k,u,A 



^ (e 'nk,.,i + e "^bk,u,2) , 
^ (e-'^6fc,.,3-e-^Hfc,.,4) . 



(139) 



(140) 



Substituting (139) and (140) into (136), we obtain 

Hj= - 2Z { 2cos ^ {bl ,^,ibk,u,i + bl ,, 3bk,u,3) + cos | {e-'hl^^^bk,u,3 + e'^4,z.,3^fc,^,i) 



k,i> 



2 COS ^ {bl^,^bk,^.,2 + bl^^^bk,^.,4) + COS ^ (e 'hl^^^bk,^.,4 + e'hl^,^bk,^.,2) 



+ 



+ isin I (e 'hl^^,bk,u,4 - e'hl^^^bk,u,i - e ''h\^,^J}k,v,3 + e'^^^s^^.^.s) |. (141) 

Here we have combined the two pairs of operators: {bk,i>,i,bk^v,3} with energies 2cos^, and 
{bk,v,2-,bk,u,4\ with energies —2 cos |. Next we observe that the diagonalization of the first two 
lines in (141) reduces to the diagonalization of only the non-diagonal parts (the second terms in 
the square brackets) which is achieved by the transformations similar to (139)-(140): 



X / • fc • fc \ 

bk,u,i = \ e~^*dk,u,i + e~'^dk,u,2) , 



\ / ■ k ■ k \ 

bk,u,3 = [e'^dk,u,i - e'^dk,u,2) , 



(142) 
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and 

1 f —i^ — ?i \ 

\ / ■ k ■ k \ 

h,u,4 = {e'^dk,u,3 - e'idk^uA) ■ (143) 
After such transformations, the Hamiltonian (141) becomes 

Hj= - Y.{ ^+4,u,ldk,u,l - e+dl^.^dk,r.^3 + i sin ^ (4,r.,l4,r.,3 " 4,z.,34,^,l) 

^-dk,u,2dk,u,2 - e-4,z.,44,:.,4 - i sin I (4,^,2^^*;,^ - 4,^,4C^fe,-^,2) 



+ 



+ 



(144) 



where 



e+ = £+(A;, I/) = 2cos - + cos -, £_ = £_(A;, z^) = 2cos — — cos — . (145) 

Thus we have obtained two independent pairs of operators: {dk^u,i, dk^u,3} with energies 

e+ and — e+, and {dk,u,2, dk,u,4} , with energies e_ and — e_. So, the diagonahzation of the 
Hamiltonian (134) is reduced to the diagonalization of two independent quadratic forms. The 
first and second lines in (144) are diagonalised respectively by the unitary transformation 



and 



dk,iy,i = cos6+Ck^i,^i-isme+Ck^t,^4, 
dk,u,3 = -i sin 9+Ck,u,i + cos 0+Ck,u,4: 



dk,u,2 = COS 9^Ck^v^2 + i sin 9-Ck,u,3, 
dk,uA = i sin 6-Cu,y,2 + cos 9-Ck,iy,3- 



(146) 



Here 9± = 9±{k, v) are determined from the relations 



tan26'± 



sm \ 



2 cos I ± cos I 



(147) 



(148) 



After this we obtain the final expression for the Hamiltonian (134) in the diagonal represen- 
tation: 



where 



Hj=Y^ [-^+A,u,lCk,v,l - S-cl,^.^Ck,u,2 + £-cl^y,zCk,u,3 + £+cl,uACk,v,A 
k,v 



£± = S±{k, u) = £± COS 29± + sin - sin 29± = 



(149) 



/ e% + sin^ 



1+4 COS^ 77 ± 4 cos - COS — . 



(150) 



Thus, the electron Hamiltonian (23) has been transformed into the diagonalised form (26) 

He = He,o + Hj= ^ Ex{k,u)cl^^^^^Ck,u,X,'r (151) 



with the energy bands (29). 
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Combining all the transformations together, we can write down the resulting unitary trans- 
formation: 

afe,i/,3 = [e-'^+Ck,u,i - e'''-Ck,y,2 - e'^-Ck,u,3) + e-'^+Ck,u,4) , 

ak,vA = ^e'^ (e^^+Cfe,^,i - e"^^" 0^,^,2 + e"^^" 0^,^,3) - e^^+ 0^,^,4) (152) 

which can be written in the general form (137). 
The inverse transformation is 

= i(e'('i'+''*)ai,M + e'('¥-'*)a,,.,, + e-<'¥-'*)a,,.,, + e-('?+''*)«.„.4 

= J(e.(^-M„,.,,,eK^-)„,..-e-.(^-).,„,3-e-'(^-M<.,„,). 

Cfc,z.,3 = ^(e*(^-^-)afc,^,i - e'( V+'^-)ofc,^,2 - e-*(^+^-)afe,^,3 + e-*(V-'^-) 0^,^,4^ , 

= I (e<^+^-)a,,,, - e<^-^+)a,,,2 + e-<^-^-)a,,,3 - e-<^+^-)a,,,4) 

(153) 

which can be written in the general form as 

Ck,u,\ = '^)ak,u,j- (154) 

13 Appendix 2. Semi- Classical Equations 

Due to the fact that Hna (55) is nondiagonal, its average value on the wavefunctions (41) 
vanishes: {^\Hna\^) = 0. The average value of the total Hamiltonian, TC = {^\H\^) = E, gives 
us the energy in the zero adiabatic approximation. The calculation of TL with (42) and (41) gives 
us the Hamiltonian functional of classical displacements and the quasiparticle wavefunction 
(Pi,j,p- So wc sec that the zero adiabatic approximation leads to the semiclassical approach which 
is often used in the description of self-trapped states. 
Calculating H, we get the Hamiltonian functional, 

n = Hph + l^ael' - J ^^*^^s{^) + xi E 

le 5 S 

+ X2M'^C^ + G2^iplip5(^)Wd^). (155) 
S 

Here Hpfi is given by (16), where the displacements and the canonically conjugate momenta are 
classical variables. The site labelling in (155) corresponds to the elementary cell with two atoms, 

based on the non-orthogonal basis vectors, and index labcles sites i,j,p- 

Prom (155) we derive the following static equations for the functions u, v, s and ip^: 

= (w + S,)^-j(^^ + ^i + ^a) + xiv{Wr + Wi + Wd) + X2VC 

+ G2['PrWr + ipiWi + ^dWd], (156) 
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I— Oi a 

V^cos{-){Wi - Wr) + cos(-)(Oi - o,) + 2^^<i 



+ kc 



sin(^)(^cos2(^)-l)(Q-a) 



+ XiiWl? - \Vr?)0^ cos(^)) + xM'' - \^r?) sm(T)(^^ cos2(t) - 1) 



,a 5 



■4' 



\/3 a 

+ — C0s(-)G2(</'Vz + V'r¥'-</cVr + </'r</^), 



(157) 



{2Wd -Wr-Wi) + ^3(0^ + 



^sin(|)(2C + a + Q) 



1 /3 



(158) 



= 
+ kc 



V3sin(^)(W^r + W^/) + sin(^)(f^^ + ili 
(2 cos(-) - - cos\-)){Cr + Ci) - cos(-)Q + ?,cos\-)C 



+ xi^sin(|)(2|^|2 + |^,|2 + |^;|2) 



+ X2((2Cos(|)--cos3(|))(|<^,|2 + |<^;|2)_cos(|)|^dp + 3cos=^(|)|^|2) 



\/3 a 

+ G2—sin(-)((^Vr + (/?;</? + </?>; + </?;>). 



(159) 



These equations were used in [23] to determine numerically the conditions for the existence 
of polaron/soliton states. 
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